Dynamical instability of a condensate induced by a rotating thermal gas 
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We study surface modes of the condensate in the presence of a rotating thermal cloud in an 
axisymmetric trap. By considering collisions that transfer atoms between the condensate and non- 
condensate, we find that modes which rotate in the same sense as the thermal cloud damp less 
strongly than counter-rotating modes. We show that above a critical angular rotation frequency, 
equivalent to the Landau stability criterion, the co-rotating mode becomes dynamically unstable, 
leading to the possibility of vortex nucleation. This kind of mechanism is proposed as a natural 
explanation f or the formation of vortices observed recently in the experiment of Haljan et al [P. C. 
Haljan et al, |cond-mat/0106362| . We also generalize our stability analysis to treat the case where 
the stationary state of the condensate already possesses a single vortex. 



The recent success at producing vortices with quan- 
tized angular momentum in trapped Bose gases g, |^ 
has raised interesting questions regarding the microscopic 
mechanisms responsible for vortex nucleation. Several 
papers [||, ||, ||, |^, ^, |^ have suggested the idea that a 
possible mechanism is the excitation of low energy sur- 
face condensate modes, which have been studied in var- 
ious experiments. This can be understood |^^, in 
terms of the well-known Landau criterion [lOl, which 
predicts the excitation of condensate modes with fre- 
quency ujm and angular momentum mh when the rota- 
tion rate fl of an external perturbation satisfies the con- 
dition ujm — mV, < 0. This defines a critical frequency, 
flc,m = oJm/m for transfer of angular momentum to the 
condensate |^, ||, where the particular mode excited de- 
pends upon the form of the "stirring" potential used. For 
example, a rotating trap anisotropy will predominantly 
couple to a quadrupole surface mode, which will become 
dynamically unstable above the Landau frequency Q,c,2- 
In related work, Sinha and Castin have found dynam- 
ical instabilities about stationary condensate solutions in 
the rotating frame jl^, and link these to vortex nucle- 
ation in a recent experiment ||l^ . Further experiments 
lend support to the role of dynamical instabilities of sur- 
face modes in the vortex formation process ^ 

In the present letter, we extend this discussion of vor- 
tex nucleation to finite temperatures and include for the 
first time the coupHng of the thermal cloud to the con- 
densate. Our calculations address the kind of experi- 
ments done at JILA in which a thermal cloud is 
set into rotational motion above Tbec by means of an 
anisotropic rotating drive. A subsequent quench through 
the Bose-Einstein condensation transition in the absence 
of the drive then leads to the formation of a condensate 
containing vortices (when the rotation frequency Vl of 
the thermal cloud exceeds a critical value). What sets 



'Present address: Electron and Optical Physics Division, National 
Institute of Standards and Technology, Geithersburg, MD 20899. 



the JILA study apart from all other recent exper- 
iments on vortices [0, ^, ^ Q is that the nucleation 
process occurs in the presence of an axially symmetric 
trap. The JILA experiment indicates that the rotating 
thermal cloud is acting as a reservoir of angular momen- 
tum that is transferred to the condensate in the process 
of vortex formation. It is clear that a breaking of rota- 
tional symmetry is required if such a transfer is to occur. 
Thermal excitation of condensate collective modes is an 
obvious candidate for this symmetry breaking. The den- 
sity fluctuations associated with such a mode are coupled 
to the thermal cloud by means of the condensate mean 
field, as well as through the collisional exchange of atoms 
between the condensate and thermal cloud. 

Of these two physical mechanisms, we focus on the lat- 
ter in this paper. The mean field coupling is also expected 
to be important, but we shall demonstrate that the col- 
lisional exchange of atoms is already sufficient to exhibit 
the dynamical instability of the condensate induced by 
the rotating thermal cloud. Generalizing recent work by 
Williams and Griflin ||l6|, we calculate the damping of 
condensate collective modes at flnite temperature due to 
the coupHng via C12 collisions between the condensate 
and noncondensate atoms. In formulating this problem, 
we assume that the thermal cloud is in rigid body ro- 
tation with angular frequency fl, as in the JILA experi- 
ments Our major result is that the damping of the 
collective modes of frequency ujm depends on the polar- 
ity m of the collective mode with respect to the rotation 
of the thermal cloud. As CI is increased, the damping of 
the TO < modes increases while that of to > modes 
decreases. At a critical angular frequency flc.m = ^m/'m 
of the thermal cloud, the to > mode becomes unstable 
in the sense that its amplitude grows exponentially with 
time. This signals a continuous transfer of angular mo- 
mentum from the thermal cloud to the condensate which 
we identify with the onset of vortex formation. 

The dynamics of the condensate wavefunction $(r, t) 
in the ZNG kinetic theory jl^ is described by a gen- 
eralized Gross-Pitaevskii (GP) equation, which at flnite 
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temperatures includes the interaction with the noncon- 
densate through a mean-field term and a non-Hermitian 
damping term due to C12 collisions that exchange atoms 
between the condensate and noncondensate. The non- 
condensate, represented by a phase-space distribution 
function /(r,p,t) obeying a semiclassical kinetic equa- 
tion is likewise coupled to the condensate through 
mean field and collisional exchange terms. The atoms are 
trapped in a static harmonic trap possessing cylindrical 
symmetry J7cxt(r) = Muj^^lp^ 4- \^z'^)/2, where the as- 
pect ratio of the trap is A = ujz/^^- The noncondensate 
atoms move in a dynamic Hartree-Fock (HF) potential 
U{r,t) — Ucyit{'r) + 2g[nc{r,t)+n{r,t)], where the interac- 
tion parameter g = iirh'^a/M, a is the s-wave scattering 
length, rtc(r, i) = |<I>(r,t)p and n{r,t) is the nonconden- 
sate density defined by n(r,i) = / c?p/(r, p, i)/(27r?i)^. 

In our idealized model of the JILA experiments, the 
collective excitations of the condensate are determined 
for an equilibrium state defined by a thermal cloud un- 
dergoing soHd body rotation at the angular frequency 17. 
For this situation, the equilibrium distribution of the ro- 
tating thermal cloud is given by 



r(r,p) 
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exp{/3[(p - Mv„o)V2M -f t/eff - Ao] } - 1 ' 

(1) 

where the noncondensate velocity field is v„o(r) — 
n X r ^ Q,p4>. The effective potential acting on the 
thermal cloud is Ucs{r) = Uo{r) - with Uo{r) 

the equilibrium HF potential. The corresponding equi- 
librium state of the condensate is given by the solution 
of the time-independent generalized GP equation. Al- 
though we are primarily interested in the initial onset 
of vortex formation, we can readily extend our analy- 
sis to a situation in which the condensate already has 
a single vortex centered on the z-axis with q units of 
quantized circulation. The equilibrium condensate wave- 
function {p, <j>, z) = V'*-''' (p, z)e"?'^ is then defined by 
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+ Y«cO + C^cxt + 5"c0 + 25^0 j = A*o'V^'^ 

(2) 

Here, the condensate velocity field is given by Vco(p) — 
{hq/Mp)4>. 

In equilibrium, the chemical potentials of the thermal 

cloud and condensate are related by /2o = Mo'"' ~ hftq. 
It is straightforward to prove that the form of /°(r, p) 
in (|l|) is a stationary solution of the ZNG kinetic equa- 
tion for the case of a static axially symmetric trap when 
the chemical potentials satisfy the above relation. Us- 
ing (|]), the thermal cloud density is given by no{p, z) — 
93/2[zo{p, z)]/ Aq, where the thermal de Broglie wave- 
length is Aq = /MkBT)^/'^ and the equilibrium 
fugacity is zo{p,z) = e^'^^'^^f'^'''^''^'^"^. As O increases, 
the aspect ratio of the noncondensate density increases 
approximately as uj^j — fi^ . 

Within our model, we neglect the dynamics of the non- 
condensate and assume the thermal cloud remains rigidly 



rotating, as described by (|^) (see also Refs. ||l6| , |l^ 
for a discussion of this type of static approximation) . In 
effect, we are ignoring the perturbation of the rotating 
normal cloud induced by the mean field of the oscillating 
condensate. Working wit h the am plitude and phase of 
the condensate <i>(r,t) = ^nc(r, t)e^^^^'*-\ the condensate 
dynamics is described by the coupled quantum hydrody- 
namic equations |^ , |l7| 



dt 
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where the local chemical potential ^c(r, t) of the conden- 



sate is 
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+ C^cxt + gnc + 25710. 



(5) 



The condensate velocity is Vc(r,t) = f?i/A/) V0(r, t). 
The source term r52(r,t) appearing in (1^) describes the 
collisional exchange of atoms between the condensate and 
thermal cloud and is given by (see also Ref. 0) 
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where the coUision time ri2 (r, t) is defined as 



^ = (2%^!' I ^Pi / / c^P3^(Pc + Pi - P2 - Pa) 



,)(l + /?)/2°/3°- 



(7) 

The local condensate momentum per atom is Pc(r,<) = 
Mvc(r,t) and the local energy is ec{r,t) — Pc{'r,t) + 
Mv^/2. The equilibrium single-particle distribution 
function /f = /"(r-jp^) is as given by (0). We note that 
the term r°2 in (g|) is the only source of damping in our 
calculations. Our static approximation for the rotating 
thermal cloud precludes Landau damping, which can be 
expected to provide additional damping of the conden- 
sate modes. 

To study the collective modes, we consider the density 
Snc{r, t) and velocity 5vc(r, t) fiuctuations of the conden- 
sate about equilibrium nc(r,t) = "-co(r) -t-(5nc(r,t) and 
Vc(r,t) = Vco(r) -H(5vc(r,t). In the Thomas- Fermi (TF) 
limit, the linearized equations of motion for the conden- 
sate fiuctuations are given by 



dSric 



+ V • {UcoSVc + ^coSuc) 



-STu 



dt 



(8) 
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It is straightforward to show that the hnearized form of 
r?2(i",t) given in (|) is 



ST, 
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dric H dvc ■ (vco - v„o) 
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(10) 



where 1/T'(r) = gUcoM /[kBTTi2{r)]. Herer°2('") is given 
by the expression in (|^), but with the dynamic or time- 
dependent condensate quantities replaced by the corre- 
sponding equihbrium ones. We remark that the second 
term in (|l^) has the form of a mutual friction term de- 
pending on the relative velocity Vco(r) — v„o(r) between 
the condensate and normal gas, which couples to the ve- 
locity fluctuation of the condensate Svc- For = and 
Vco — 0, the results in reduce to those in Ref. |Q. 

In the TF approximation, the condensate density is 



given by rico (r) = {h^q^ -Ue^t- "^gno - 2 - 
plify the calculation, we neglect the effect of the noncon- 
densate mean field on the condensate profile |jl^. In ad- 
dition, we neglect the effect of the vortex core on the over- 
all shape of the condensate. This approximation is valid 
when ^^/Rq ^ 1, where ^ is the healing length and Rq 
is the Thomas-Fermi condensate radius ||2^. We there- 
fore approximate the condensate density as n^o — (ptf — 
Ucxt)/9, where f iTFjT) = huj_L[15XNc{T)a/dA_]'^^^ /2 and 

In an axially symmetric trap, the density and ve- 
locity fluctuations of the collective modes take the 
form Snc{r,t) = z)e*(™'^-'^'"*) and Svc{r,t) = 

(?i/M)V(56'™(/9,z)e*(™'^-'^™*). Substituting these expres- 
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(11) 



(12) 



where S ee g[\/ ■ {nco'\7) -m^rico/ p'^]/M and V = fiq/Mp. 
These equations can be combined into a single equation 
for 5nm, following the procedure used in Ref. We 
relate 59m to Srim, making use of the following approxi- 
mation to ( |l2|) 
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This is vaHd for p :s> consistent with neglecting the 
effect of the vortex core on the static condensate density 
proflle. We then obtain a single wave equation for the 
condensate density fluctuation for mode m 



{uim + S + Sq + iSr)Snr, 
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~ _ 2mhq 

« = UJrnMp^ 



2gn 



cO 



Mp 



d_ 

dp 



(14) 



(15) 



and 



St = (^rn - mil) — . 
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In order to obtain an approximate solution of (|lj) , we 
treat the effect of the vortex Sq and the damping Sr 
perturbatively ||l^ by expanding in the zeroth-order TF 
solutions given by S6ns — —ojgSns ||2^. Eq. ( |l4| ) can 
then be simplifled to | [l^ 



(17) 



where the spatial average of an operator x is defined by 
(x) = / drx(r)(5n|(r)/ J dr(5n|(r). The frequency yi?^ is 
the collective mode frequency of a vortex state 



(18) 



where lus 
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is the TF frequency for g = 0, with 
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"^^'p^ M p3 \dp p 



_ (19) 
The solution of (|lj) to lowest order in the damping is 

LUm — i^rn — i^rn , where the damping rate is given by 



r(.)-i/l\('i_!!^ 

" 2\r' V ,„(9) 



(20) 



We note that the finite T result in (Q) is formally iden- 
tical to the T = result in Ref. |^^. The only difference 
is the use of Nc{T) for the number of condensate atoms. 
As discussed in Refs. ||2^, collective modes with op- 
posite polarity m are split in frequency since the vortex 
circulation breaks the azimuthal symmetry. Modes rotat- 
ing in the same sense as the vortex (to > 0) are shifted 
higher in frequency and modes rotating against the fiow 
of the vortex (to < 0) are shifted down, in agreement 
with experiments H, |l^ . 

The expression for the damping (|20| ) of the collective 
modes of a condensate interacting with a rigidly rotating 
thermal gas is our main new result. Interestingly, it pre- 
dicts that modes of opposite polarity to are affected by 
the rotating thermal cloud in quite distinct ways. The 
damping of modes that rotate against the flow of the 
thermal cloud (to < 0) increases linearly with angular 
frequency Q, while the damping of modes that rotate in 
the same sense as the thermal cloud (to > 0) decreases 
with increasing fl. This prediction should be easily ob- 
servable. For TO > excitations, the damping in ( [20| ) 
goes to zero when il reaches a critical value defined by 



,(?) 



(21) 



For r2 > rii'm, the damping changes sign, indicating the 
onset of a dynamical instability (i.e., Sum ~ e"'"''"™!*). 
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The result in ( ^l| ) is equivalent to the usual Landau cri- 
terion for the excitation of surface modes by an external 
anisotropic perturbation as discussed in Ref. [Q. In con- 
trast, our finite-temperature mechanism is effective even 
in an axisymmetric trap. Physically, one would expect 
all of the surface modes to be thermally occupied, so that 
for a noncondensate rotation rate of fl, instabilities can 
be induced in any of these modes so long as > fli'^m- 
In other words, at least one mode will be unstable when 
fl exceeds a critical frequency defined by the minimum 
of a plot of iJ"m I'm versus m [Q. This critical rotation 
frequency is generally larger than that required to ensure 
thermodynamic stability of a vortex state js) , so that it is 
likely that the primary role of surface mode instabilities 
is to facilitate tunneling of one or more vortex lines into 
the condensate bulk. This has recently |Q been demon- 
strated using energy arguments for the particular case of 
a m = 2 quadrupole mode. When a single vortex 
is already present, our result (18) shows that the crit- 
ical rotation frequencies increase due to corresponding 
upward shifts in the to > frequencies |2^. Again, 
taking the to = 2 mode as a concrete example, our TF 
analysis predicts = uj^/\/2 without a vortex, while 

^c^2 = t^_Lv/(l + 2A2)/2 for a g = 1 vortex state. After 
nucleation, the vortices will eventually equilibrate into 
a thermodynamically stable lattice configuration, as dis- 
cussed in im,!!. 

We emphasize that the present paper only evaluates 
the damping from C12 collisions in the "collisionless" or 



mean field region. We expect that Landau damping, 
which would arise if we included thermal cloud fiuctu- 
ations |6|, would also exhibit the same dynamical 

instability for > Jli^m- The damping term in ( |lo| ) can 
be interpreted as a kind of mutual friction, although it is 
quite different from the usual mutual friction arising in 
the two-fiuid hydrodynamic domain, such as discussed in 
Ref. 

In summary, stimulated by recent work at T = on 
the nucleation of vortices by a rotating anisotropic har- 
monic potential, we have generalized the discussion to 
finite temperatures and considered the case where the 
driving field is a rigidly rotating thermal cloud . Our 
calculation is explicitly based on the collisional exchange 
coupling between this rotating thermal cloud and the 
condensate which gives rise to mode damping. This 
damping can change sign when the thermal cloud angu- 
lar frequency Vl reaches a critical value given by the Lan- 
dau criterion for excitation of collective modes of angular 
momentum m along the z axis. Importantly, our result 
also applies to the stability of a singly quantized vortex 
against excitations of surface collective modes, where it 
is associated with dynamic nucleation of additional vor- 
tices. 
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a critical reading of an early draft of this paper. This 
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of Japan. 
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